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The steady-state nucleation rate and flux of composite nucleus at the saddle point is studied by extending 
the theory of binary nucleation. The Fokker-Planck equation that describes the nucleation flux is derived 
using the Master equation for the growth of the composite nucleus, which consists of the core of the final 
stable phase surrounded by a wetting layer of the intermediate metastable phase nucleated from a metastablc 
parent phase recently evaluated by the author [J. Chem. Phys. 134, 164508 (2011)]. The Fokker-Planck 
equation is similar to that used in the theory of binary nucleation, but the non-diagonal elements exist in 
the reaction rate matrix. First, the general solution for the steady-state nucleation rate and the direction of 
nucleation flux is derived. Next, this information is then used to study the nucleation of composite nucleus 
at the saddle point. The dependence of steady-state nucleation rate as well as the direction of nucleation 
flux on the reaction rate in addition to the free-energy surface is studied using a model free-energy surface. 
The direction of nucleation current deviates from the steepest-descent direction of the free-energy surface. 
The results show the importance of two reaction rate constants: one from the metastable environment to the 
intermediate metastable phase and the other from the metastable intermediate phase to the stable new phase. 
On the other hand, the gradient of the potential <£> or the Kramers crossover function (the commitment or 
splitting probability) is relatively insensitive to reaction rates or free-energy surface. 
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I. INTRODUCTION 

Nucleation is a very basic phenomena which plays a 
vital role in various material processing in industry rang- 
ing from steel production to food and beverage indus- 
tries 1 . Recently, researchers have focused on the nu- 
cleation of complex materials which are relevant to our 
daily life 1 . The nucleation of such complex materials can 
also be complex and may involve intermediate metastablc 
phases 2-5 . These metastable intermediate phases indi- 
cate a composite nucleus that consists of a core of stable 
new phase surrounded by a wetting layer 3 of intermedi- 
ate phase, which is nucleated from the metastable parent 
phase. 

The author 6 recently examined the thermodynamics 
and free-energy surface of a critical nucleus using cap- 
illarity approximation when an intermediate metastablc 
phase exists The author discovered the possibility of such 
a composite critical nucleus at the saddle point of the 
free-energy surface 6 . In order to study the free-energy 
surface of nucleation, we utilized a three-phase system 
and assumed the existence of an intermediate metastable 
phase in addition to an unstable parent phase and a sta- 
ble new phase. We chose a composite binary nucleus that 
consists of the stable new phase and the metastable inter- 
mediate phase. This binary nucleus is completely phase 
separated, and the metastable phase is segregated on to 
the surface. By studying the free-energy surface of the 
two-component system in the two-dimensional space , we 
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can easily visualize the free-energy surface and are able 
to locate the saddle point that corresponds to the criti- 
cal nucleus. Even though such a composite nucleus has 
long been predicted theoretically 7-9 , it has been proved 
explicitly that the composite nucleus corresponds to the 
saddle point only through the use of a computer simu- 
lation 10 . Our calculations 6 clearly demonstrate that by 
using an analytical model the composite nucleus really 
corresponds to the saddle point of the free-energy sur- 
face. A similar composite nucleus at the saddle point of 
the free-energy surface is considered in the problem of 
deliquescence 11-13 . 

From the free-energy surface, the author could eas- 
ily calculate the so-called minimum-free-energy path 
(MFEP) 14 . However, it has been recognized for more 
than thirty years that the steepest-descent direction on 
the free-energy surface 15 that corresponds to the MFEP 
does not necessarily indicate the real nucleation pathway 
in binary nucleation 16-18 . This is because the asymme- 
try of the reaction rate will deflect the direction of nucle- 
ation flux from the steepest-descent direction. Likewise, 
the free-energy surface alone may not determine the nu- 
cleation pathway for the composite nucleus. The rea- 
son for this is that the asymmetry between the reaction 
rate from the metastable parent phase to the intermedi- 
ate metastable phase and that from the intermediate to 
stable phase will play some role to determine not only 
the nucleation rate but the direction of nucleation flux. 
Therefore, the MFEP determined by the steepest-descent 
direction of free-energy surface 14 will not necessarily cor- 
respond to the real nucleation pathway. 

This study focuses on the steady-state nucleation rate 
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and flux at the saddle point for the composite nu- 
cleus through the theory of binary nucleation 19,20 . We 
must first determine the Fokker-Planck or the Zeldovich- 
Frenkel equation 1 ' 21-23 from the Master equation for the 
composite nucleus (Section II) . The theory of nucleation 
rate and flux for the binary nucleation when the reac- 
tion rate matrix is non-diagonal is re-constructed in Sec- 
tion III and in the Appendix A and B as the theory of 
Trinkaus 19 is too condensed to study various aspects of 
nucleation of multicomponent systems. In section IV the 
theory developed in Sections II and III will be applied to 
the composite nucleus. An example to illustrate the pro- 
cess will be provided in order to show how the reaction 
rate will influence the nucleation rate and the direction 
of nucleation flux at the saddle point for the compos- 
ite nucleus. Section V will contain the conclusion of the 
study. 



II. FOKKER-PLANCK EQUATION FOR THE 
COMPOSITE NUCLEUS 



In order to study the nucleation kinetics of the com- 
posite nucleus, the model shown in Fig. 1 is considered. 
The model consists of a core of the stable new phase 
(number of molecules ni) surrounded by an intermedi- 
ate metastable phase (number of molecules n 2 ) nucle- 
ated in the metastable parent phase 3,6,9,10 . Nucleation 
rates k + and k~ are the reaction rates between the sta- 
ble new phase and the metastable intermediate phase. 
Nucleation rates a + and oT are the reaction rate or at- 
tachment rate from the metastable parent phase to the 
intermediate metastable phase. 

A similar model has been used to study the parti- 
tioning transformation 24-26 where the coupling of the in- 
terfacial and long-range diffusion fluxes is considered to 
study the time-dependent precipitate nucleation 27 , void 
formation in irradiated metals 28 , and nanocrystal forma- 
tion from metallic glasses 29 . In these linked- flux mod- 
els 24-26 , in contrast to our model 6 , the metastable inter- 
mediate phase does not really represent the thermody- 
namics phase. Rather, it models a transition zone due 
to the diffusion of monomers 24,25 . Therefore, the nucle- 
ation rate a + does not represents the attachment rate of 
monomer but is simply a parameter which approximately 
represents diffusion. In our model, however, the interme- 
diate phase represents the metastable thermodynamics 
phase and a + is the attachment rate of the monomer. 

The Master equation for the time-dependence of 
the number of clusters / (n-i, n 2 , t) that consists of ri\ 
molecules of the stable phase in the core and n 2 molecules 
of the intermediate metastable phase in the surrounding 
wetting layer is written generally 1,15 ' 21 in the form 




Metastable intermediate phase 



Metastable parent phase 



FIG. 1. A composite critical nucleus model that consists of 
a stable new phase (number of molecules ni) surrounded by 
an intermediate metastable phase (number of molecules 712) 
nucleated in the metastable parent phase 6 . Nucleation rates 
k + and k~ are the reaction rates between the stable new 
phase and the metastable intermediate phase. Nucleation 
rates q + and a~ are the reaction rate or attachment rate from 
the metastable parent phase to the intermediate metastable 
phase. 



where 



J ni = « + {ni,n 2 ) f (ni,n 2 ,t) 

- k~ (ni + 1, n 2 — 1) / (?ii + 1, n 2 - 1, t) , 

Jni-l = « + ("1 -1,^2 + 1)/ (™1 - 1,«2 + 1,*) 

- k~ (ni,n 2 ) f(n 1) n 2) t) , (2) 
J n2 = a + (m, n 2 ) / (ni, n 2 , t) 

~ a~ (m, n 2 + 1) / (ni, n 2 + 1, t) , 
J n2 -i = a + (ni,n 2 - 1) f (ni,n 2 - l,t) 

- a~ (ni,n 2 )/(ni,n 2 ,t). 



Using the detailed balance condition, 



K (ni + 1, n 2 - 1) = k + (m, n 2 ) 



/ eq (ni,n 2 ) 



/eq ("1 + l,n 2 - 1) : 



- 1 1 -n + 1 \ /eq (ni,n 2 ) 

a (ni,n 2 + l) = Q! + (ni,n 2 )-— ? — , (3) 

/eq (ni,n 2 + 1) 



where the equilibrium cluster distribution / oq (n) is given 
by the usual Boltzmann distribution 



/eq (n) = fo exp (-/3G (n)) , 



(4) 



and G (n) is the work of cluster formation for a clus- 
ter with composition n = (ni,n 2 ) and /3 is the inverse 
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temperature. Then, Eq. (1) can be written as 

d/Qi,n 2 ,£) _ 
dt 

- k + (ni,n 2 ) / eq (ni,n 2 ) 

7 (n 1 ,n 2 ,t) _ /(ni + l,n 2 - l,t) 
/ cq (ni,n 2 ) / eq (ni + l,n 2 - 1) 
+ k+ (m - l,n 2 + f ) / eq (ni - l,n 2 + L 
" / (rti - l,n 2 + M) _ / (ni,n 2 ,t) 
./oq(«l - 1,"2 + 1) feq(ni,n 2 ) 

- a+ (ni,n 2 )/ eq (ni,n 2 ) 
f{ni,n 2 ,t) /(ni,n 2 + l,t) 



./ eq (ni,n 2 ) / e q(ni,n 2 + l) 
+ a+ (ni,n 2 - 1) / eq (ni,n 2 - 1) 
~f{ni,n 2 -1,*) f(n 1 ,n 2 ,t) 



Jeq(ni,n 2 -1) f eci (ni,n 2 ) 
In the continuum limit, we have 



g/(n) 

at 



<9J ni 9J, 



+ 



9ni <9n 2 



-divJ, 



(5) 



(6) 



where the components of the nucleation flux J are given 

by 



L ani t)n 2 
J -* + f l — \ 

_ K+/ e 1-^ + ^ 

eq \ 9ni <9n 2 



and 



$ (ni,n 2 ,f) 



f(n 1 ,n 2 ,t) 
/ eq (ni,n 2 ) ' 



(7) 



(8) 



(9) 



and simplified notations k+ = k+ (m,n 2 ), a + = 
a+ (ni,n 2 ) and / eq = / eq (ni,n 2 ) and $ = <&(m,n 2 ,t) 
are used. Equation (6) is called the Fokker-Planck 30 or 
the Zeldovich-Frenkel equation 22,23 . 

Equations (7) and (8) can be put in the form of a 
matrix equation 



/eq (n)fl(n)V„$(n), 



(10) 



where .7 and V„ are row vectors, and R is a symmetric 
square matrix defined through Eqs. (7) and (8). The 
two fluxes J ni and J„ 2 are linked by the non-diagonal 
rate matrix R. Eq. (10) is originally derived by Russell 24 
though the analysis of the nucleation flux at the saddle 
point was mostly qualitative and numerical and did not 
present various analytical formulas which will be derived 
in this paper. 

A general theory of barrier crossing by the decay of 
metastable state has been formulated by Langcr 31 using 
the Fokker-Planck equation more than forty years ago 



by generalizing the Kramers' theory 32 of barrier crossing, 
and has been extended by several workers 33,3 to study 
the reaction and nucleation pathways in various problems 
not limited to binary nucleation. This Kramers-Langer- 
Berezhkovskii-Szabo (KLBS) theory 31-33 considers a gen- 
eral problem of escape probability from a potential well 
of the metastable state. In this KLBS theory, the nu- 
cleation flux J is the probability flux and $ is called 
the Kramers crossover function 34 or the commitment or 
splitting probability 35 which represents the probability 
to cross the barrier from the given point in the phase 
space. 

In this paper, we adopt the theory of binary nucle- 
ation 15-21 since it does not need to characterize the ini- 
tial metastable state by a potential well. Also, the 
connections to the original Master equation are more 
transparent since the elementary process of nucleation 
is a discrete process of attachment and detachment of 
monomers. It is well recognized 16 that we must return 
to the original discrete Master equation from the Fokker- 
Planck equation when we study the nucleation pathway 
of post-critical nucleus. In any case, the formula based 
on KLBS theory 31-33 reduces to that used in the theory 
of binary nucleation 19 when the thermal distribution in 
the metastable potential well is incorporated 36 . We will 
extend the theory 15-21 developed for binary nucleation 
to study the nucleation of composite nucleus in the next 
section. 



III. MAGNITUDE AND DIRECTION OF STEADY 
STATE NUCLEATION FLUX AT THE SADDLE POINT 

From this point, we will concentrate on the steady 
state nucleation flux at the saddle point n* so that the 
nucleation flux is a time independent constant vector 
(div J = 0) and the rate matrix R = R (n* ) is also a 
constant matrix at the saddle point. 

The general theory of nucleation flux of binary nucle- 
ation will be reformulated when the rate matrix R is 
non-diagonal so that we can study the nucleation of com- 
posite nucleus. Most of the previous authors who studied 
the nucleation flux of binary nucleation have studied the 
cases where the rate matrix is diagonal 16,20 . To exemplify 
this, we introduce the rotation matrix of the coordinate 
given by 

(W? -sin0\ 
v ' \sm6 cos 9 J v ' 

which is the orthogonal transformation in vector space, 
and study the steady-state solution df /dt = of the 
Fokker-Planck equation (6) given by 



^ n J — Qri\ Jni "F d n2 J n 



0. 



(12) 



where we have used a column vector for the differen- 
tial operator div and the superscript T which means the 
transpose. 
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Using the rotation matrix in Eq. (11), we can intro- 
duce a new coordinate e v whose origin is displaced at 
the saddle point n* and rotate the unit vectors e ni and 
e„ 2 by an angle 9 through 

e n = V(9)e v , or e v = V T (6) e n . (13) 

Then the displacement (column) vector from the saddle 
point n* 



**-»- »•-(£:) 



(14) 



will be transformed into the displacement n through 

An = V(9)n, or n = V T (9) An. (15) 

The direction of nucleation flux given by Eq. (10) can be 
studied by inspecting the appropriate rotation of coordi- 
nate. 

A. Direction of steepest-descent of free-energy surface 

The steepest-descent direction of the free-energy sur- 
face, which is the direction of minimum free-energy path 
(MFEP) 14 , at the saddle point can be studied by ex- 
panding the work of cluster formation G (n*) around the 
saddle point n* as 15 



G (n) = G* + ^An T GAn 



(16) 



using the displacement (column) vector defined in 
Eq. (14), where 



with 



G = 



Gi 



Gn G\2 
Goi G 



22 



d 2 G 
drudrij ' 



(17) 



(18) 



is the Hessian matrix of the free energy at the saddle 
point n* , and G* = G (n*). By using the rotation matrix 
V (9) defined in Eq. (11), we can diagonalize the matrix 
G and obtain the negative eigenvalue that corresponds 
to the steepest-descent direction as 

Ai = ((Gii + G 22 ) - V / (G 11 +G 22 ) 2 -4dctG^ /2, 

(19) 



B. Direction of gradient of $ 

Next, we will consider the direction of the gradient of 
$, which has been studied intensively by Wilemski 20 and 
others 37 ' 38 . This function $ is also known as the Kramers 
crossover function 34 or the commitment or splitting prob- 
ability 35 in the general theory of barrier crossing, where 
$ plays the role of the probability to cross the barrier. 
We will extend the work of Wilemski 20 since we are deal- 
ing with nucleation by an attachment and detachment 
process. Now, we rotate (Am, An 2 ) axis by lo and intro- 
duce a new coordinate (Ci, C2) through 



C = V T (to) A n . 



(21) 



Then, the component and Jj 2 are related to J ni and 
J„ 2 through 



J C = V T (u)J n . 



Similarly we obtain 



V c = V T (lo) V T 



(22) 



(23) 



which is the simplest transformation of covariant vec- 
tor 30 ' 39,40 when the coordinate transformation is given 
by a linear orthogonal transformation Eq. (11) as the 
gradient of a scalar is a covariant vector. The covariant 
formulation of the Fokkcr-Planck equation 30 ' 39 will be 
useful when we consider more general coordinate trans- 
formations 40 ' 41 . 

The nucleation flux in Eq. (10) is now transformed into 

= -/eq («) V T (to) RV (w) V c $ (24) 

at the saddle point. 

By calculating the matrix products V T KV, Eq. (24) 
is written explicitly by 

J Cl = -/ eq (n) S ( W ) || + / eq (n) S(w)J^, (25) 
J C2 = f cq (n) H M || - /«, (n) T ( w ) (26) 
where 

S (u) — Ru cos 2 bj + R 2 2 sin 2 lo + 2i? i2 cos uj sin u, 



(27) 



where detG = GnG 22 - G( 2 < at the saddle point. 5 ( w ) = (jj n _ ij 22 ) coscj sinw + R 12 (sin 2 cj - cos 2 w) 



The rotation angle is chosen such that Ai becomes 
negative and A 2 becomes positive. Then the rotation 
angle 9 is given by 



T (u>) = Ru sin 2 uj + i? 22 cos 2 lo — 2i? i2 cos u sin to 



, /, (G 22 - Gn) , 
tan ^ 2G 12 ± 




(28) 



(29) 



+ 1, (20) 



where + sign must be chosen when Gi 2 < otherwise 
— sign must be chosen. The steepest-descent direction 
is the direction of rji axis defined through from Eq. (13) 
with 9 given by Eq. (20). This is the direction of the 
eigenvector of negative eigenvalue Ai of the matrix G. 



are the matrix elements of V T RV. 

We then choose the new coordinate £ such that c\ 2 $ = 
0. After manipulating Eqs. (25) and (26) using the 
same procedure originally developed by Langer 31 , which 
is briefly summarized in Appendix A, we find 



tan to = s ± \/ s 1 + r, 



(30) 
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where 



R22G 



s = 



22^22 



RuG 



11 



2 (-R12G11 + R 22 G\ 2 ) 
R11G12 + R21G22 
R12G11 + R22G12 ' 



(31) 
(32) 



are slightly different in definition from those used by oth- 
ers 16 ' 20 since we have off-diagonal elements R12 — R 2 \- 

Eq. (30) reduces to Eq. (20) that determines the 
steepest-descent direction of free-energy surface when 
^11 = -R22 and i?i2 = -R21 = 0. Again, + sign must 
be chosen when R\ 2 Gn + R22G12 < otherwise — sign 
must be chosen. Equation (30) with + sign also reduces 
to the well-known result 16,20 when i?i 2 = R21 = 0. The 
+ sign was chosen by previous authors because they 16 ' 20 
considered the case when G\ 2 < 0. 

Equation (A3) can now be written as 



d 2 <S> 



-CiA 



with 



A = L (w) /E (w) . 
Since the boundary condition is given by 

$(Ci) = 1, as Ci -> -00, 
$(Ci) = 0, as Ci -> 00, 
the solution of Eq. (33) is given by 

*(&) = ^erfc 




(33) 

(34) 

(35) 
(36) 

(37) 



where erfc is the complementary error function. There- 
fore, the gradient of $ is given by 

d<5> 



with 



cos oje ni + sinwe^ 



(38) 



(39) 



from Eq. (21) and the rotation angle u is given by 
Eq. (30). It is easily confirmed by direct calculation that 
the direction of given by Eqs. (30)-(32) is in fact 
the direction of the eigenvector of the matrix GR (Ap- 
pendix) . 

Equation (37) indicates that the steady state cluster 
distribution / (n\, n 2 ) is smaller than the equilibrium dis- 
tribution / eq (711,712) since <f> < 1 and / = / oq /2 at the 
saddle point. Also, $ = /// eq becomes a universal func- 
tion along (1 axis. Numerical results of vapor to liquid 
binary nucleation 37 ' 38 also suggest that whole contour of 
<f> = /// eq in two-dimensional (ni,n 2 ) space can also be 
relatively insensitive to the materials considered. There- 
fore the universality of Eq. (37) could be used as a guide 
to deduce the reaction coordinate on the free-energy sur- 
face and the generalized coordinates that correspond to 
n\ and n 2 using various techniques of computer simula- 
tions 42 - 45 . 



C. Direction of nucleation flux J 

Thirdly, we will derive the direction of the nucleation 
flux J following the work of Wilemski 20 . Now, we rotate 
(Arii,A7i2) axis by <j> and introduce a new coordinate 
(fii,fi 2 ). Then, the nucleation flux given by Eq. (10) will 
be transformed into 



J„=/eq (») V T (0) RV (</>) 



(40) 



similar to Eq. (24). 

Using a procedure similar to that used in the previous 
subsection and Appendix A 20 ' 31 , we obtain the relation- 
ship between </> and to given by 



tan < 



R12 + R22 tancj 
i?n + i?i2tanw' 



(41) 



which reduces to the formula derived previously 16 ' 20 
when i?i2 = R 2 i =0. 

Similarly, the nucleation flux J is given by 



(42) 



with 



J Ml =/ eq (n) (detfl) [a^)]- 1 ^ 



exp - 



Ad 2 , , 



(43) 



where a (</>) is given by 

a ((f)) = [(R 2 n sin 2 (f> + R% 2 cos 2 <$) 

-2R 12 (Ru + R 22 ) cos (/.sin + R 2 12 ] 1/2 .(44) 
and its direction is 



'Ml 



cos <pe n , + sin <pe ri 



(45) 



with the rotation angle <j> given by Eq. (41). 

Equation (41) tells us that the direction = 
(cos 4>, sin of the nucleation flux e m (cx e^) and the 
direction vj x — (cos cj, sin lj ) of the gradient of $ (e^) 
is related through cx cx Rv < ^ 1 , which again con- 
firms the statement in the Appendix that the direction 
of the eigenvector of the negative eigenvalue Ai of the 
matrix GR represents the direction of the gradient of $ 
and that of the matrix RG represents the direction of 
the nucleation flux J at the saddle point. 



D. Nucleation rate 

Finally, we will consider the nucleation rate. Since 
the general theory was already given by Trinkaus 19 , 
we will simply apply his theory to the two-dimensional 
case in this subsection. In Appendix B, we give 
a detailed derivation of the formula for the nucle- 
ation rate by Trinkaus 19 and the relationship between 
the Trinkaus' theory and Kramers-Langer-Berezhkovskii- 
Szabo (KLBS) theory 31 - 33 . 
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To simplify the notation, we will introduce matrices 

r = R^ 2 = ( ri1 n A (46) 
\r 2 i r 22 J y ' 

and 



r =^ = (Kl£) (17) 

defined in Eq. (BIO). We also introduce the rotation 
of coordinate V (a) in Eq. (fl) with rotation angle a 
that corresponds to the orthonormal transformation in 
Eq. (B14). 

Now, the nucleation flux is given as Eqs. (42) -(45) is 
also written by Eqs. (B18) and (B19). Therefore 



and 



Ad 2 = 



/3Ai fa(<j>) 



(48) 



(49) 



since 



and 



(50) 



detR \p{4>), 
\e tl \=&etR 1 ' 2 / P { ( l)) 
P(4>) = [On sin 2 4> + r 2 22 cos 2 0) 

1/2 

-2ri2 (rn + r 22 ) cos sin ^ + r 2 2 ] , (51) 

using the direction <j> of the nucleation flux. 

Using the same procedure as that used in section III A, 
we have one negative (Ai) and one positive (A 2 ) eigen- 
values 



Ai = ^(r u + r 22 ) - y/(T u + r 22 ) - 4detrj /2 < o, 

(52) 

a 2 = ((rn + r 22 ) + yjpu + r 22 ) 2 - 4dctr^) /2 > o, 

(53) 

and the rotation angle a is given by 



r 22 — Tn 

tan a = — — ± 



22 — 1 n 



2r 



12 



2T 



12 



+ 1 



(54) 



where + sign must be chosen if T 12 < otherwise — sign 
must be chosen. Since, the direction of nucleation flux 
in Eq. (B19) must be parallel to e Ml in Eq. (45), we 
have 



tan 6 — 



T12 + ?"22 tan a 
rn + r\2 tana 



(55) 



which will be reduced to the form obtained by Wilemski 20 
as rn = R\{ 2 < 
(ri2 = r 2 i = 0). 



1/2 1/2 

as r n = i? n and r 22 = R 22 when R 12 = i? 21 = 



Finally, the nucleation rate is given by 



foe 



-PG\ 



I detR 



Ail 



(56) 



from Eq. (B21). 

In the next section, we will study the three angles 6, 
uj and <f> which represent the steepest-descent direction 
on the free-energy surface, gradient of $ (commitment 
probability), and the nucleation flux J. These directions 
are characterized by 9, uj and <f> which correspond to the 
directions of the eigenvector of the negative eigenvalue of 
the matrix G, GR and RG, respectively 31,33 . 



IV. APPLICATION TO THE COMPOSITE NUCLEUS 

In this section, we will apply the theory of nucleation 
flux of binary nucleation at the saddle point developed 
in sec. Ill to the problem of nucleation flux of composite 
nucleus considered in sec. II. In this case, the reaction 
rate matrix is non-diagonal, 



R = 



c+ + a+ 



(57) 



as given in Eq. (10). Similar models with non-diagonal 
reaction matrix were considered by Russel 24 as a linked 
flux model and by Djikaev 11 to study the nucleation 
in deliquescence using the theory of Melikhov ct al. 46 . 
However the latter author used the theory 46 which re- 
duces the two-dimensional problem to an approximate 
one-dimensional problem. Furthermore, only a numerical 
result of deliquescence path is given and a detailed the- 
oretical analysis is missing 11 ' 24 . Here, we will study the 
two-dimensional problem directly using rigorously formu- 
lated theory without resorting to an approximation. 

Since, the eigenvalues of T are the same as those of 
the product RG and GR, they can easily be obtain from 
Eqs. (52) and (53) explicitly as 

Ai= \ [(a+G 22 + K + (Gn + G 22 - 2G 12 )) 



(a+G 22 + k+ (Gn + G 22 - 2G i2 )) 2 - 4a+fc+dctG 



(58) 



A 2 = - [(a+G 2 



«+ (Gi 



G 2 



2G 12 )) 



+ \j {a+G 22 + n+ (Gn + G 22 - 2G i2 )) 2 - 4a+fc+detG 



(59) 



where Ai < and A 2 > since detG < at the saddle 
point. Using these eigenvalues Ai and A 2 , the nucleation 
rate is given by 



/ = foe 



|Ai| 



(60) 



Nucleation rate and flux of composite nucleus 



7 




FIG. 2. Relation between uj and <f> for various relative reaction 
rate a = a + /k + . The angle <f> is an increasing function of w 
starting from <j> = —45° at uj = 0°. The straight line indicates 
the line when <f> — u. Below this line, the direction of the 
nucleation current is deflected from the direction of grad$ to 
the right. 



from Eq. (56) since deti? = a + K + . The nucleation rate 
depends not only on the rate a + from the metastable 
parent phase to the metastable intermediate phase but 
also on the rate k + from the metastable intermediate 
phase to the stable new phase. Therefore, even when 
a + 3> k + , the nucleation problem does not reduce to one- 
component unary problem 19,36 ' 46,47 . Since there is only 
one saddle point, a formula for the successive nucleation 
rate 6,48 ' 49 cannot apply. 

The steepest-descent direction 8 of free-energy sur- 
face is given by Eq. (20), which is determined solely by 
the free-energy surface near the critical point given by 
Eq. (16). The direction uj of the gradient of $, however, 
is given by a more complex formula in Eq. (30) because 
the effect of the anisotropy of reaction rate cannot be ig- 
nored. Using the explicit form of reaction rate in Eq. (57) 
we have the direction d> 



tan i 



-k + + (k + + a + ) tanw 
«+ — k + tanw 



(61) 



once we know uj from Eqs. (30). 

Figure 2 shows the relationship between uj and <f> for 
various relative reaction rate a = a + /n + . The angle 
<j) is an increasing function of the angle uj starting from 
<j> = -45° at uj = 0° from Eq, (61). The direction <j> of 
nucleation current J changes from right {<p < uj) to left 
(<j) > uj) of the direction uj of the gradient of <E> as the 
angle uj is increased. As the angle uj is further increased, 
the angle <j> becomes <f> = 90° when uj = 45° from Eq. (61) 
and continue to increase. Therefore, as the direction uj of 
the gradient $ increases, the nucleation flux at the saddle 
point turns toward the 712-axis and the component n 2 of 
the intermediate metastable phase rather than n\ of the 
stable new phase tends to grow. 



When a < 1 (a + < k + ), the incoming rate a + from the 
metastable parent phase to the intermediate metastable 
phase (712) is slower than the internal transformation rate 
k + from the intermediate metastable phase to the new 
stable phase. Then the stable new phase (ni) tends 
to grow faster than the intermediate metastable phase 
(712) j and the direction <f> of nucleation flux remains to- 
ward the direction of n\ axis and the inequality <j> < uj 
follows as long as uj is not as large as 45° (see Fig. 2, 
a = 0.1,0.7). When a > 1, the incoming rate a + is 
faster than the internal transformation rate k + and the 
intermediate metastable phase grows faster than the new 
stable phase. Then the direction (f> soon becomes larger 
than uj as uj is increased from (see Fig. 2, a = 2.0, 5.0). 

In order to study the relationship and the absolute 
magnitude of various angles 9, uj, and 0, we have to spec- 
ify the free-energy surface G (n) . To enable the qualita- 
tive idea of the composite nucleus, we will give an exam- 
ple to illustrate the concept and study a model system 
specified by the free-energy landscape 

G(n) = G* 

+ \g 22 



1 + ^§12 ) An i + A ™2 + 2gi 2 AniAn 2 



(62) 



with gi2 < 0. In this case G\\ = G22 (l + 2/v / 3gi2) , 
G12 = G22gi2, and we find from Eq. (20): 



tan 6=—=, or 9 = 30°. 
V3 



(63) 



In this model, only gi 2 < is the parameter that charac- 
terizes the free-energy landscape. In order to make the 
origin (Ani,An 2 ) = (0,0) the saddle point, detG must 
be negative and gi2 must satisfy gi2 < — 1/V3 = —0.577. 

In Figs. 3 and 4, we show contour plots of typical free- 
energy landscape of our model in Eq. (62) when gi 2 = 
—0.7 and gi2 = —5.0. The saddle point is located at 
the origin (0, 0). Although the direction of the steepest- 
descent is fixed to 30°, the saddle point and the valley 
that lead to the saddle becomes wider and shallower as 
the parameter g'12 becomes more negative (cf. Figs. 3 
and 4). 

Figure 5 shows the angles uj and <f> as the functions 
of the relative reaction rate a. Since we always have 
RvtGu + i? 22 Gi2 = -K+G u + (k+ + a+) G 12 < be- 
cause G11 > and G12 < 0, we took + sign in Eq. (30). 
We can easily show from Eqs. (30), (61) that 



tancj 



1 - 



G22 — G 



12 



G11 + G2 



2Gi 



as a 



Therefore 



uj — > 45°, as a + 



and 



tan < 



G11 — Gi 



G 2 



Gn 



as a 



•0+ 

(64) 



(65) 



(66) 
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-1 -0.5 0.5 1 



FIG. 3. An example of the model free-energy landscape of 
composite nucleus when 6 = 30° defined by Eq. (62) when 
gi2 = —0.7. The origin is the saddle point shifted to (0,0). 
The steepest-descent direction is always fixed to 6 — 30° in- 
dicated by a pink arrow from the saddle point indicated by a 
pink point. The direction ui of the gradient of $ in Fig. 4 is in- 
dicated by a red thin (a = 0.0) and a thick (a — 7.0) arrows. 
The direction <j> of the nucleation flux J is not shown as it 
stays near the steepest-descent direction 9 = 30° represented 
by a pink arrow (see also Fig. 5). 




-1 -0.5 0.5 1 



FIG. 4. The same as Fig. 3 when gi2 = —5.0. Now the value of 
gi2 becomes more negative than that in Fig. 3 and the saddle 
point and the valley that leads to the saddle become wider. 
The direction ui of the gradient of $ in Fig. 5 is indicated by 
green thin (a = 0.0) and thick (a — 7.0) arrows of solid line. 
The direction <j> of nucleation flux J is indicated by green thin 
(a = 0.0) and thick (a — 7.0) arrows of broken line. 



70 
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a 

FIG. 5. The angle uj and (j> as the functions of the relative 
reaction rate & for various free-energy parameters gi2 when 
61 = 30°. 



from Eq. (61). The angle uj that specifies the direction 
of the gradient of $ always starts from uj — 45° at a — 
irrespective of the shape of the free-energy surface G (n) . 

As the relative reaction rate a is increased, the incom- 
ing flux from the metastable parent phase to the sur- 
rounding wetting layer of the intermediate metastable 
phase increases. Then the intermediate metastable com- 
ponent ri2 of the composite nucleus increases and the 
nucleation flux tends toward direction of rii axis and the 
angle </> increases. In Fig. 5, the angle u, which is the 
direction of gradient of $, starts from u> = 45° and ro- 
tates clockwise towards n\ axis, while the angle cA, which 
is the direction of nucleation flux, starts from the right- 
hand side of the direction 9 and rotates anti-clockwise 
and turns to the left-hand side of the direction of 9 as 
the relative reaction rate a is increased. 

In Figs. 3 and 4, we show the direction of w and <p at 
a = 0.0 and 7.0 by thin and thick arrows using the same 
color as those used in Fig. 5. In Fig. 4, the arrow of solid 
line which represents the direction of u> of the gradient of 
$ rotates clockwise, while that of broken line which rep- 
resents the direction of </> of the nucleation flux J rotates 
anti-clockwise as a is increased. It can be shown analyt- 
ically that both the directions uj and <f> become the same 
as the steepest-descent direction 9 = (f> = uj = 30° at 
a = 2/V3 (see Fig. 5). Whengi 2 = — 1/V3, the direction 
cj) is fixed to the steepest-descent direction <f> = 6 = 30° 
because the free-energy landscape becomes the corridor 
of valley rather than the saddle. 

In Fig. 6, we show the direction u> and <p as the function 
of the relative reaction rate a when the steepest-descent 
direction is 9 = 60° using the model free-energy land- 
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scape 

G{n) = G* 
1 



gi2 Anf + An" 



2gi 2 AmAn 2 
(67) 



As the relative reaction rate a is increased, the angle 
u) starts from to = 45° and rotate clockwise towards 
ri\ axis again. The angle oj is relatively insensitive to 
the parameter gi 2 of the free-energy surface not only in 
Fig. 6 but also in Fig 5. Several authors observed similar 
behaviors of oj and $ in vapor-liquid nucleation 20 ' 37 ' 50 
which are relatively insensitive to the system consid- 
ered. In our case, the decreasing oj indicates that the 
size distribution function / will be characterized mainly 
by the number of molecules ni in the final stable phase 
(/ (711,112) ~ / as a is increased. The size distribu- 
tion / (ni, 71.2) is flat along n% axis for fixed n\ and many 
clusters with the same n\ and different ni coexist. 

On the other hand, the angle <j> starts from the angle 
(j> > 60° calculated from Eq. (66) that is larger than the 
steepest-descent direction 9 = 60° and continues to in- 
crease. When gi2 = — \/3, the direction <f> is fixed again 
to the steepest-descent direction <j> = = 60° because 
the free-energy landscape becomes corridor of a valley. 
In Fig. 6 the angle oj and (j> will never coincide. As ex- 
pected, the direction <f> of the nucleation flux is closer 
than the direction oj of the gradient of $ to the direc- 
tion 9 of the steepest-descent of the free-energy surface. 
Our results in Figs. 5 and 6 clearly show that the dif- 
ference between the steepest-descent direction 9 and the 
nucleation flux direction cf> depends not only on the dif- 
ference of the reaction rates represented by a but also 
to the shape of the free-energy surface as pointed out by 
Wyslouzil and Wilemski 51 . 




FIG. 6. The same as Fig 5 when 9 = 60°. 

The deviation of the direction of the nucleation flux 
from that of the steepest-descent direction of free-energy 
surface has been pointed out repeatedly in different con- 



text from condensed matter 17-19,24 to liquid-vapor sys- 
tems 15 ' 16 ' 20 . However those authors except Trinkaus 19 
used various approximations, in particular, the reaction 
rate matrix is assumed to be diagonal. Therefore, their 
conclusion remained qualitative. In contrast, we have 
used the theory 31-33 which includes the non-diagonal ele- 
ment of the reaction rate matrix R to study the direction 
of nucleation flux. Therefore, our result will be useful to 
study the nucleation flux of more complex situation and 
will contribute to deeper understanding of a more com- 
plex nucleation process in the future. 



V. CONCLUDING REMARK 

In this paper, we have studied the nucleation flux of 
composite nucleus that consists of a stable new phase sur- 
rounded by a wetting layer of intermediate metastablc 
phase nucleated from the metastable parents phase. 
By using the results of exact formulation of nucleation 
flux 31-33 , we have shown that the nucleation rate is de- 
termined not only from the reaction rate from the inter- 
mediate metastable phase to the stable new phase but 
also from that from the metastable parent phase to the 
intermediate metastable phase. We have also shown that 
the direction of nucleation flux deviates significantly from 
that of the steepest-descent direction of the free-energy 
surface. 

Kinetics of composite nucleus has been considered by 
several authors in different contexts. Kelton and cowork- 
ers 25 ' 26 have extended the theory of linked flux theory 
of nucleation by Russell 24 and studied the partitioning 
transition where the interplay of diffusional and interfa- 
cial processes are linked. In this case, the nucleus consists 
of a core of new stable phase surrounded by a layer of pro- 
nucleus zone. They 25,26 did not show any analysis of the 
nucleation flux since they do not appear to be interested 
in the size distribution of nucleus which is directly com- 
parable to a real experiment. In fact, Kelton 2 applied 
the theory to the oxygen precipitation in silicon and ob- 
tained a precipitation density that agrees fairly well with 
their experiment. Russell 24 showed only numerically that 
the nucleation flux will deviate from the steepest-descent 
direction of free-energy surface. Very recent study by Pe- 
ters 52 , however, does not use the wetting layer as a model 
of diffusion but considers the coupling of diffusion and 
nucleation directly. His results clearly demonstrated that 
the nucleation pathway is influenced by the diffusion rep- 
resented by along-wavelength fluctuation and also by the 
initial state and questioned the adequacy of the compos- 
ite nucleus model. Recent studies 53 ' 54 using the dynami- 
cal density functional theory also pointed out the impor- 
tance of long-wavelength fluctuation as the pre-nucleus 
state and questioned the composite nucleus model as a 
model of two-step nucleation. 

A similar composite nucleus that appears in the prob- 
lem of deliquescence has been considered by Djikaev 11 . 
In this model 11 the wetting layer is not merely a dense 
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pre-nucleus zone 2 but is rather a true thermodynamic 
liquid phase. Therefore his model is closer to ours 6 since 
we consider the wetting layer as a metastable thermo- 
dynamic phase. For a realistic model, he 11 showed only 
numerically that the deliquescence path, which is the di- 
rection of nucleation flux, deviates from the direction of 
the equilibrium path of the steepest-descent on the free- 
energy surface. In contrast, our analytical analysis is 
not restricted to specific materials and will be useful for 
future extension and applications to various complex ma- 
terials. 

It must be noted, however, our analysis is completely 
confined to the steady-state process. Therefore, transient 
properties 25 ' 26,36,47 and the problem of time scale 11,46 
are out of scope for this study, which can only be ex- 
amined numerically by solving coupled Master equa- 
tion 25,26,55 . Also, our analysis assumed that the nucle- 
ation flux goes through the saddle point. Saddle point 
avoidance 19,36,51,56,57 will be important if the anisotropy 
of the reaction matrix R is large or the ridge between sad- 
dle point is low, which can occur at high temperatures 
or near the spinodal point. In such a case, ridge-crossing 
rather than saddle-crossing may occur. Then the nucle- 
ation flux will spread over the whole phase space and the 
picture used in this study may break down. In such a 
case, the growth of composite nucleus will be expected 
to be more complex. 

Finally, our analysis will be useful in the study of nu- 
cleation by attachment and detachment of molecules as 
well as the study of various barrier crossing phenomena 
such as chemical reactions and conformation transforma- 
tion of molecules. Newly developed various computer 
simulation techniques make it possible not only to ex- 
plore the free-energy landscape 42,44,45 but also to extracts 
the reaction rate matrix 34,58 and identify reaction coor- 
dinate 34, 59-62 . Therefore, our analysis will be useful to 
deduce important parameters from those numerical data 
extracted from simulations. 

Our analysis will also be useful even when multiple 
saddles exist 6, 63-66 and multiple pathways of nucleation 
coexist. In this case, the nucleation flux may fork into 
multiple streams and an analogy to the electric current in 
a circuit with parallel and series connections of conduc- 
tors or resistors will be useful 67 as long as the fluxes go 
through those saddles. Our results will be valid if some 
of those multiple saddles correspond to the composite 
nucleus. 
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Appendix A: Determination of the direction of gradient of 

$ 

In this appendix, we will determine the £ coordinate 
system such that the direction of £i axis corresponds 
to the direction of gradient of $ using the same pro- 
cedure originally developed by Langer 31 and used later 
by Wilemski 20 . Now, we will choose the new coordinate 
C such that = 0. Then, 

J Cl =-/«,£ (w)0 Cl $, (Al) 
J C2 =/ eq A(a;)%$, (A2) 

and the steady-state condition = can be written 

as a differential equation 

rf 2 $ d<E> 

where 

P((1) . SM (^)^ AM (^) ( , (A4) 

Then, using the transformation (21), 

In jcq = j3 ^G* + ^C T V T (w) GV (w) Cj + constant, 

(A5) 

and after calculating the matrix product V T GV, and the 
derivatives 9 In f eq /d(i and 9 In f eq /d(2, Eq. (A4) can be 
written as 

p(Ci)=Cii(w) + C2M(w), (A6) 

where 

L (u>) = -/3E (u>) (Gn cos 2 u + G 22 sin 2 ui + G i2 sin 2w) 
+/3S (w) ((G 22 — Gn) sin w cos co + G i2 cos2u;) , 

(A7) 

M (u>) = -/3E (u>) ((G 22 - Gn) sinwcosw + Gi 2 cos2cj) 
(w) (G 22 cos 2 cj + Gn sin 2 co — Gi 2 sin 2ui) . 

(A8) 

Since it is assumed that $ does not depend on £ 2 , £ 2 - 
dependence of p(Ci) in Eq. (A6) must be eliminated. 
Therefore, we have to set M = in Eq. (A8). Using 
Eqs. (27) and (28) and setting M = in Eq. (A8), we 
obtain the rotation angle co that satisfies Eqs. (30)-(32). 

Appendix B: Saddle point nucleation of multicomponent 
systems 

In this appendix, we provide a summary of the theory 
of nucleation rate for general multicomponent systems 



Nucleation rate and flux of composite nucleus 



11 



by Trinkaus 19 since his presentation is too condensed to 
indicate the detail of the derivation of various formulas. 

The starting point is the the steady-state nucleation 
flux represented by a column vector J defined by 



J=-R(n)f cq (n)V 
in Eq. (10). Using the unit vector 



f(n) 

/eq (n) 



(Bl) 



(B2) 



\ e n„ 



along the orthogonal axis (m, n 2 , . . . , n m ), we can write 
the nucleation current explicitly in the vector notation 

J = J 1 e ni +J 2 e, l2 ^ h J m e, lm = J T e n = e^J. (B3) 

The equilibrium cluster distribution / eq (n) is given by 
the usual Boltzmann distribution similar to Eq. (4). 
Therefore, the flux is written as 

J = R (n) cxp (-0G (n)) V n exp (+0G (n)) / (n) , 

(B4) 

which will be solved with the boundary conditions 

(B5) 
(B6) 



/(|n|^0)^/ eq (|n|^0), 
/ (|n| -> oo) -> 0. 



In order to study the nucleation flux in Eq. (B4) near 
the saddle point, we expand free energy G (n) around the 
saddle point n* according to Eq. (16). Next, in order to 
eliminate the anisotropy introduced by the reaction rate 
matrix R, we introducing a new variable u through a 
linear transformation 



R}/ 2 



An, 



(B7) 



where R}l 2 is the symmetric square root matrix of R. 
Since R}l 2 is not an orthogonal transformation, the base 
vectors are transformed from e„ to e v according to 



R}l\ 



(B8) 



It must be noted that the new basis e v are neither or- 
thogonal nor normalized. 

In the new coordinate system i/, the steady-state nu- 
cleation flux in Eq. (B3) is written as 



J = -e 



(G*+» T r„/2) v ^ e +0(G*+„ T rv/2) f ^ 



from Eq. (B4), where 

T = R}' 2 GR 1 ' 2 . 



(B9) 



(BIO) 



Since the matrix T defined in Eq. (BIO) is a symmet- 
ric matrix because G and R 1 / 2 are symmetric, we can 
diagonalizc T using an orthogonal matrix V as 



V T TV = A, 



(Bll) 



where A is a diagonal matrix of the eigenvalues Ai given 

by 



Tvi = XiVi 



l,...,n, 



(B12) 



where Vi is the normalized (column) eigenvectors. We 
choose Ai < as a single negative eigenvalues. Other 
eigenvalues A2 , . . . , A m are positive since the point n* is 
the saddle point and dctG < 0. The orthogonal (or- 
thonormal) matrix V is given by 



V = (vi,v 2 , ...,«„ 



(B13) 



using the (normalized) eigenvectors v^. 

Using the this orthogonal transformation V from v to 
i by 



I = V T u 1 

which corresponds to the rotation of axis 
e € = V T e„, 



(B14) 



(B15) 



we can diagonalize matrix T. Then, the steady-state nu- 
cleation flux is given by 

J = _ e -/3(G*+rA«/2) v ^(G'+eM/2)f (|) ) 

(B16) 

which must be a time-independent constant. 

It is reasonable to assume that the nucleation flux 
goes to the direction that corresponds to the neg- 
ative eigenvalue Ai of T in Eq. (B12). Then, the multi- 
dimensional problem in Eq. (B16) reduces to the one di- 
mensional problem 



-,3(G*+£ T AS/2)_^_ e/ 3(G*+£ T AS/2) j ^ 



(B17) 



which can be integrated by £1. Using the boundary con- 
dition Eqs. (B5) and (B6) and the expansion (16) around 
the saddle point, we arrive at the equation 



^1 = /eq O) 



2tt 



cxp 



?\£ 2 
2 lCl 



and 



(Bis) 



(B19) 



These results were originally derived by Trinkaus 19 . We 
should note that the basis vector is not normalized 
(|e^J 7^ 1). Therefore its magnitude contribute to the 
net magnitude of nucleation flux J in Eq. (B19). 

The nucleation rate / can be calculated by integrating 
the nucleation flux inm-1 dimensional hypersurface 
in n, which can be transformed into the integration by 
£2, • • • , im'- 

/oo />oo 
■ ■ ■ / Jfc {£) detR 1/2 d&d& ■ ■ ■ dtm, (B20) 
-OO J — OO 
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where detii 1//2 is the Jacobian to preserve the volume 
element. Using Eq. (B18) in Eq. (B20), wc find 



I= PM detR V* foe -?G*^ 



n 



2tt 



(B21) 



Since 



[J A 4 = dctA = dety T rF = detF T i? 1/2 Gil 1/ V 

i=l 

= (detR^^detG, (B22) 



we can rewrite Eq. (B21) in the form 



7 = 



/3|Ai|/2tt 
v/dct [I3G/2tt 



■foe 



-pa* 



(B23) 



derived originally by Langer 31 and re-derived later by 
Trinkaus 19 . 

Finally, we will note the connection to the more general 
KLBS theory 31-33 . It is easy to show from Eqs. (BIO) and 

(B12) that the direction = vje u = (R}-I 2 v\) e n of 
the nucleation flux is the direction of the eigenvector of 
the negative eigenvalue Ai of the matrix RG. Therefore, 
the eigenvector f + defined by Berezhkovskii and Szabo 33 
represents the direction of nucleation flux. 

Similarly, it is easy to show from Eqs. (BIO) and 
(B12) that the direction of the eigenvector of the negative 
eigenvalue Ai of the matrix GR represents the direction 
of the gradient of the stochastic separatrix defined by 
$ = i/2 33 ' 35 - 58 . Therefore, the eigenvector e+ defined 
by Berezhkovskii and Szabo 33 represents the direction 
of V$ at the saddle point in Eqs. (38) and (39) in 
our nucleation problem. 
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